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ABSTRACT 

We employ an all-sky map of the anomalous microwave emission (AME) produced by component 
separation of the microwave sky to study correlations between the AME and Galactic dust properties. 
We find that while the AME is highly correlated with all tracers of dust emission, the best predictor 
of the AME strength is the dust radiance. Fluctuations in the AME intensity per dust radiance are 
uncorrelated with fluctuations in the emission from polycyclic aromatic hydrocarbons (PAHs), casting 
doubt on the association between AME and PAHs. The PAH abundance is strongly correlated with 
the dust optical depth and dust radiance, consistent with PAH destruction in low density regions. 
We find that the AME intensity increases with increasing radiation field strength, at variance with 
predictions from the spinning dust hypothesis. Finally, the temperature-dependence of the AME per 
dust radiance disfavors the interpretation of the AME as thermal emission. A reconsideration of other 
AME carriers, such as ultrasmall silicates, and other emission mechanisms, such as magnetic dipole 
emission, is warranted. 


1. INTRODUCTION 

High sensitivity, full-sky observations of the far- 
infrared and microwave sky from WMAP and Planck 
have pushed studies of the Cosmic Microwave Back¬ 
ground (CMB) to a regime in which contamination from 
Galactic foregrounds has become a key uncertainty in 
the analysis. Understanding the physical nature of the 
foreground components and producing better models of 
each are essential for mitigating this uncertainty. 

The anomalous microwave emission (AME) is perhaps 
the least well-understood of the foreground components. 
Discovered as a dust-correlated emission excess peak- 


ing near 30 GHz (Kogut et al. 

1996 de Oliveira-Costa 

et al. 

1997 Leitch et al. 1997|, 

A1V1E is often ascribed 


small dust 

grains (Draine & Lazarian 

1998a|b; Hoang 

et al. 2010 

Y sard V Verstraete 

2010 

; Sif 

sbee et al.|201l|, 


sion comp onent peak s between ~ 20 and 50 GHz (e.g 


Planck Collaboration et alj2014b) and has an e missivity 
peFTT^r^33TlT^^iysF^^n?II _1 at 30GHz dDobler 

™ ' '■ it 


et al. 2009 
tion et al. 


Tibbs et al 
2014b| 


c n. 


2011| Planck Collabora" 


,,20101 

Polycyclic aromatic hydrocarbons 
(PAHs), which give rise to prominent emission features in 
the infrared, are considered natural carriers of the AME 


due to their small size and apparent abundance ([Draine 
& Lazarian||1998a). 
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Theoretical spinning dust SED templates based on a 


brandon.s.hensley@jpl.nasa.gov 


PAH size distribution that reproduces the infrared emis¬ 
sion features have been successful in fi tting observati ons 
of the AME both in the Galaxy flMiville-Deschenes et al 


2008 


Hoang et al.||20l"H |Planck Co llabora tion et a T. 


2015bi and in the sole extragalacti c AME detecti on 


the star-forming £ 

galaxy NGC6946 (Murphy et al. 2010 

Scaife et al. 
The AML 

2010 

J has 

Hensley et al. 
been observec 

2015JF - 

to correlate well with 


the PAH emission features in the infrared. lYsard et al.l 


(2010) found that, over the full sky, the AME was more 
correlated with emission at 12 /im than with 100 /im. 
Likewise, the AME in the dark cloud LDN 1622 was bet¬ 
ter correlated with the 12 and 25 pm emission than with 
either the 60 or 100 /im emission (Casassus et al. 2006). 

^ A i\ /fur' ^ ^^ i^ iu i '5 i^ ^ 


Howev er, studying the AME in a sample of 98 Galactic 


clouds, Planck Collaboration et al. (2014b) found no sig¬ 
nificantmHerencesHDetweerTTEeTz, 2^7^b 5 and 100 /am 
emission in their correlation with the AME. Analysis of 
both the Perseus molecular cloud (Tibbs et al.|2011) and 
the HII region RCW175 found no compelling link be¬ 
tween the PAH abundance and the AME. Likewise, the 
link between the AME and PAH abundance determined 
from dust model fitting has proven tenuous in NGC 6946 
dHensley et al.||2015|. 

In addition to spinning ultrasmall grains, grains con¬ 
taining ferro- or ferrimagnetic materials are also pre¬ 
dicted to radiate stro ngly in the microwave and may con¬ 
tribute to the AME ([ Draine & Lazarian 1999; Draine & 
Hensley 2013). Draine V Lazarian (1999) argued that the 


spinning dust and magnetic emission mechanisms could 
be distinguished by observing the AME in dense regions 
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where PAHs would likely be depleted due to coagulation. 
If the AME per dust mass is constant across both dense 
and diffuse regions, then spinning dust emission would 
be disfavored. 

In this work, we test the spinning PAH hypothesis 
using new full- sky observatio ns of t he infrared and mi¬ 


crowave sky. Planck Collaboration et al. (2015b) decom¬ 
posed the Planck: sky into foreground "components mak¬ 
ing use of both the 9-year WMAP data and the Haslam 
408 MHz survey. The combination of these data allowed 
the low frequency foreground components - primarily 
synchrotron, free-free, and AME - to be disentangled, 
producing a full-sky map of the AME. 

All-sky WISE observations at 12 fi m, a tracer of PAH 
emission, have natural synergy with the AME map, al¬ 
lowing us to test at high significance the link between 


also calculate the uncertainty in this quantity assuming 
Gaussian errors and ignoring the uncertainty on the peak 
frequency of both components. 

In addition to the AME map, we employ the param¬ 
eters from the thermal dust fit to compute the 353 GHz 
dust optical depth r 353 . Although this parameter was 
also an explicit data product^ of_the fu ll-sky modified 
blackbody fits performed by [Planck Collaboration et al. 


(|2014a|), we prefer to us e the results from |Planck Col- 


laboration et al. (2015b) for several reasons. Eirst, the 


fits were performed on more data and with more de¬ 
tailed treatment of calibration and bandpass uncertain¬ 
ties. Second, the fits were performed at 1 ° versus 5' 
resolution, mitigating the effects of cosmic infrared back¬ 
ground (CIB) anisotropies. Finally, a single-temperature 
modified blackbody model has been demonstrated to be 


ing by Planck Collaboration et al. (2014a) permits deeper 

HFI frequencies to 100 /am ( 

Meisner & Finkbeiner 2015 

exploration into the dust and environmental parameters 

Planck Collaboration et al. 

2015b). Thus, we prefer to 


This paper is organized as follows: in Section[ 2 ]we sum¬ 
marize the data sets used in our analysis; in Section [ 3 ] we 
describe the correlations predicted by our current under¬ 
standing of the spinning PAH hypothesis; in Section [4] 
we present the relationships between environmental and 
dust properties and AME as derived from the data; in 
Section [5] we discuss the implications of these results on 
spinning dust theory in particular and AME modeling in 
general; and we summarize our principal conclusions in 
Section [ 6 [ 

2. DATA 

2 . 1 . Planck Foreground Separation Maps 
Com bining full-missio n all-s ky Planck observations 


(Planck Collabo ratio n et al. 
W1V1AP data (Benn ett et al. 
408 MH z su rvey 


2015a) with the 9-yr. 
2013) 


and the Haslam 


iHasIam et al.| |l!J82|), |Planck Col- 
2015b)' used data from 32 different 


laboration et al. 
detectors spanning a range of 408 MHz to 857 GHz 
in frequency to perform foreground component separa- 
tion within the Bayes i an Commande r analysis framework 


(Eriksen et al. 2004 


2006, 2008). They constructed 


theoretically-motivated models lor the frequency depen¬ 
dence of each component - including the CMB, syn¬ 
chrotron, free-free, thermal dust, and AME - while mini¬ 
mizing the number of free parameters needed and simul¬ 
taneously fitting for calibration offsets. Using a Gibbs 
sampling algorithm, they determined the best-fit values 
for each model parameter on a pixel-by-pixel basis. 

In this work, we focus primarily upon the resulting map 
of the AME. This component was modeled by the sum 
of two spinning dust spectra wit h fixed spectral shape as 


determined by the S pDust code (Ah-Ha’imoud et al. 2009 
Silsbee et al. 2011 ), but differing amplitudes and peak 
frequencies. One of the spectra was required to have a 
spatially fixed peak frequency, fit to be 33.35 GHz, while 
the other peak frequency was allowed to freely vary from 
pixel to pixel. Thus, the data products consist of an 
amplitude for each AME component, the peak frequency 
of the spatially varying component, and the uncertainty 
of each for every pixel on the sky. 

To facilitate comparisons with the literature, we cal¬ 
culate the sum of the two components at 30 GHz. We 


of the dust emission spectrum only. We find 7353 derived 
in this way to be on average 10 % lo wer t han that reported 


by Planck Collaboration et al. (2014a). 

Additionally we employ the parameter maps that char¬ 
acterize the free-free, synchrotron, and CO emission. As 
free-free and CO emission cannot be fit reliably in regions 
of low surface brightness, we exclude all pixels with emis¬ 
sion measure of 0.1 cm -6 pc or less and with no fit CO 
emission when analyzing the former and latter, respec¬ 
tively. 

The Commander parameter maps have a resolution of 
1° (F WHM) and are pixellated with HEALPix (Gorski 
et al.||2005|) resolution of N sid e = 256. For our analysis, 
we downgrade these maps to 7V S ide = 128, corresponding 
to pixels of about 27’.5 on a side. 

2 . 2 . Dust Parameter Maps from Modified Blackbody 

Fitting 

A key limitatio n of the t hermal dust fits performed by 
Planck Collaboration et al. (2015b) is the omission of 
any data at higher frequency than the Planck 857 GHz 
band. Without information on the Wien side of the dust 
emission spectrum, it is difficult to constrain the dust 
temperature and luminosity. 


Thus, we emplo y the full-s ky parameter maps of|Planck 
Collaboration et al. (2014a), who fit a modified black¬ 
body model of the"dust emission to the 2013 Planck 857, 
545, and 353 GHz data as well as 100 /am data from IRAS. 
For the latter , t hey e mploy both t he rep rocessed IRIS 
map (Miville-Deschenes & Lagache 2005) as well as the 
map of Schlegel et al. (1998). The fits were performed 
at 5' resolution and yielded an estimate of 7353 , the dust 
radiance IZ = f I^ ust diy, the dust temperature Ik, and 
the dust spectral index (3 = d lnr/d \nv for each pixel. 
We employ 7X^, and /3 from these fits in our analysis, 
but use t 353 from the Commander fits as discussed above. 

We use the parameter map pixellated with HEALPix 
resolution of N S id e = 2048, then smooth with a Gaus¬ 
sian beam of FWHM 1° and downgrade the resolution 
to N side = 128. 

2.3. WISE 12pim Map 

WISE observed the full sky in f our in frare d bands 
- 3.4, 4.6, 12 , and 22 /am (Wright et al. [2010). The 
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12 /a m channel (W 3), sensitive from 7 to 17 /im, captures 
the strongest of the infrared emission features associated 
with PAHs and thus traces the population of small dust 
grains. The full-sky WISE 12 pm imagin g has been re¬ 
processed by Meisner & Finkbeiner ([2014) to isolate dif¬ 
fuse emission from Galactic dust at 15” resolution. 

Because interplanetary dust models are insufficiently 
accur ate to rem ove zo diacal light from the W 3 data, 
Meisner & Finkbeiner (2014) adopted a zero level based 
on the Planck 857 GHz map, essentially replacing 12 /im 
modes on scales larger than 2 ° with a rescaling of 
857 GHz dust emission. Since we wish to study real fluc¬ 
tuations in PAH emission per FIR dust radiance on 1° 
scales, this zero level procedure has the potential to ar¬ 
tificially suppress part of our signal. 

Our analysis relies on quantifying the PAH abundance 
in a given pixel through the ratio of the W3 emission to 
the FIR du st ra diance 7£, a quantity we designate /pah 
(see Section |3. 2 1) . Co mputing /pah via simple division of 
the Meisner Finkbeiner (2014) W 3 map by 7 Z is prob- 
lematic because this method is sensitive to the Planck 
857 GHz-based 12 /im large-scale zero level. 

We therefore employ an alternative approach to com¬ 
pute /pah? leveraging the excellent angular resolution of 
WISE and Planck to render the large-scale zero level 
of the W3 map irrelevant. Our principal analysis em¬ 
ploys Aside = 128 HealPix pixels of ^0.5° on a side. To 
determine /pah, we subdivide each of these pixels into 
approximately 120 subpixels of 2.5' on a side and then 
derive the best-fit correlation slope of W3 emission ver¬ 
sus 7 Z across the 120 subpixels. To infer the optimal 
/pah value within a single HEALPix pixel, we employ the 
model 


{yiv 


,12 a 


= ^i x /pah + C 


(i) 


12 

where, in the ith subpixel, {yl u ) i is the W3 map 
smoothed to 5' resolution and high-pass filtered at 15 / , 
and IZj is the Pla nck radiance map (Planck Collabora¬ 
tion et al. 2014a) high-pass filtered at 15'. The high- 


pass filtering serves to remove large-scale modes from 
the comparison, and the offset C is a nuisance parame¬ 
ter. Because linear regression can become sensitive to the 
choice of fitting methodology in the limit of poor signal- 
to-noise (see, e.g.,|Hogg et al.|2010), we restrict analyses 
of these /pah correlation slopes to HEALPix pixels with 
very strong linear correlation between (z//^) 12 ^ m and 7£, 
using a threshold of Pearson r > 0.6. This threshold 
roughly corresponds to a 10 % fractional uncertainty on 
the inferred /pah- 

The /pah map derived in this way is qualitatively simi¬ 
lar to that obtained by simply dividing the 12 /am map by 
7but displays a factor of ^2 broader range of fpAH val¬ 
ues. Our basic conclusions are unaffected by the choice 
of /pah map, and we employ the map obtained with 
Equation [l] as the default. 

Only in the case of Figure [7] do we require a map of 
the 12 /im emission itself, rather than /pah- For this 
purpose, we generat ed a custom version of the jMeisner 
& Finkbeiner (2014) W3 map by replacing only modes 
larger than 4.5° with a rescaling based on the Planck 
857 GHz emission. This better preserves real fluctuations 
in 12 /im emission per IZ on ~ 1 ° scales. 

Both our final maps of WISE 12 /am and /pah are pixel- 


lated with HEALPix resolution of N s id e = 128. 

2.4. IRAS 100/am Map 


follow Planck Collaboration et al. 

(2014a] 

) in employing 

the IRAS 100 /am map of Schlegel 

et al. ( 

1998). We do 


Des chen es fc Lag ac he 2005) o n sma ll scales as was done in 


Planck Collaboration et al. (2014a). We note that using 
the IRIS map instead in our analysis somewhat degrades 
the tightness of the correlation between the 100 /am emis¬ 
sion and the AME, particularly in regions of low surface 
brightness. The effect is small and has limited impact on 
our conclusions. 

We smooth the IRAS 100 /am map of [Schlegel et al.| 


(1998) to a FWHM of 1° and pixellate with N S id e = 128. 

2.5. Masks 

Zodiacal light in the ecliptic plane can dominate the 
signal from PAH emission at 12 /am even on small an¬ 
gular scales, thereby biasing our primary PAH tracer. 
We therefore exclude all pixels within five degrees of the 
ecliptic plane (8.7% of the sky) to mitigate this effect. 

Likewise, artifacts from moon contamination are also 
present in the 12 /am map. We therefo re mask all pixels 
flagge d for moon contamination by [Meisner & Fink beiner 
(2014), totaling 16% of the sky. Since our analysis re- 


quires that both IZ and 7/ 2Axm measure true dust emis 
sion, and since these quantities will closely correlate 
when they do, we require every pixel in our analysis to 
have a Pearson r > 0.6 between IZ and 7 / 2 /xm across 
its 2.5’ subpixels as described in Section |2.3j This cut 
alone eliminates 42% of the sky, mostly at low surface 
brightness and high Galactic latitude where WISE lacks 
sufficient sensitivity to measure diffuse dust emission. 

Though nominally a full-sky mission, IRAS did not 
have 100% sky coverage . We therefore mask th e regions 
IRAS did not observe (|Wheelock et al.||1994|), totaling 
2 % of the sky, as the 100 /am emission in thes e regions 


esti mated using lower resolution DIRBE data (Schlegel 
et al.|[l 998) may not be reliable 


To mitigate the effects of point sources, we emplo' 
the Planck point sour ce ma sks in intensity at ea ch HF] 
and LFfl band (|Planck Collaboration et al.| [2015af7 
These masks eliminate point sources with signal-to-noise 
greater than 5 at resolution N s [de = 2048. We downgrade 
these masks to N s id e = 128 by rejecting any pixel con¬ 
taining a point source, resulting in 39% of the sky being 
masked due to point sources. 

Finally, we use the Planck Galactic plane mask cov¬ 
ering 1 % of the s ky based on the 353 GHz HFI dats[^] 
(Planck Collaboration et al. 2015a). This eliminates the 
regions in the Galactic plane with the most complicated 
and highest intensity emission and where the relatively 
simple foreground models are most likely to break down. 

After applying the masks described above, 51,579 pix¬ 
els covering 26% of the sky remain. The total mask is 

1 http://irsa.ipac.caltech.edu/data/Planck/release_2/ 

ancillary-data/masks/HFI_Mask_PointSrc_2048_R2.00.fits 
^ http://irsa.ipac.caltech.edu/data/Planck/release_2/ 

ancillary-data/masks/LFI_Mask_PointSrc_2048_R2.00.fits 
° http://irsa.ipac.caltech.edu/data/Planck/release_2/ 
ancillary-data/masks/HFI_Mask_GalPlane-apo0_2048_R2.00. 
fits 
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Fig. 1.— Top: The Mollweide projection full-sky map of 1Z as derived by |Planck Collaboratio n et al.| ( |2014a| > with the mask described in 


Section 
same mask 


2.5 


overlaid in gray. The unmasked area comprises 26% of the sky. Bottom: The full-sky map of /pah = (vlv) 12 /7l with the 

overlaid. 
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illustrated in Figure [l] on the full-sky TZ map. This com¬ 
bination of masks is used in all analysis, although we 
discuss the sensitivity of our results to various additi onal 


masking (e.g. masking the Galactic plane) in Section 4.6 


3. TESTS OF AME THEORY 
3.1. AME Theory 

The AME has been suggested to be electric dipole 
emission from ultrasmall, rapidly-rotating grains that 
have been torqued up t hrough interactions with both the 
gas and radiation field (Draine fc Lazarian 1998b; Hoang 


et al.|2010||Ysard & Verstraete 2010; Silsbee et al.|201ip. 
The observed peak frequency of the emission requires 
that the grains be small (radius a < 10 A), leading to a 
natural association with the PAHs that produce emission 
features in the infrared. 

If so, we might expect a linear relation between the 
total PAH surface density Epah and the AME intensity. 
Adopting an empirical 30 GHz spinning dust emissivity 
of 3 x 10 18 Jysr 1 cm 2 H 1 , and assuming this corre¬ 
sponds to typical Galactic value s of M^/M# = 0.01 and 
Epah/E^ = 0.046 ( Draine fc Li||2007| ), we would expect 


rAME _ p o 

V30GHz — u -° 


-'PAH 


M© kpc 


-2 


Jy sr 


( 2 ) 


In the context of the spinning dust model, environmen¬ 
tal factors can influence both the peak frequency of the 
emission and the emissivity itself. The importance of col¬ 
lisions with ions depends on the fractional ionization and 
the charge state of the ultrasmall grains. In regions with 
very intense radiation fields, drag forces from IR photon 
emission become important. Damping by the rotational 
emission generally causes the ultrasmall grains to have 
sub-LTE rotation rates. 

The electric dipole moment distribution of the dust 
population will also strongly influence the emissivity, 
though we have no a priori estimates of the systematic 
variations of this quantity from one region to another. 

We expect spinning dust emission to be relatively in¬ 
sensitive to the strength of the radiation field, which 
is an important source of excitation only in fairly ex¬ 
treme environments such as reflection nebulae and photo¬ 
dissociation regions (PDRs). However, the emissivity 
per unit gas mass should increase with increasing local 
gas density, which may correlate with the radiation field 
strength. 

The impact of these fa ctors on the 30 GHz AME fl ux 
density was estimated by Draine & Lazarian (1998b) to 
be only about a factor of two between the Cold Neu- 
tral Medium, Warm Neutral Medium, and Warm Ion¬ 
ized Medium environments. In our study, these effects 
are mitigated by the low angular resolution of the maps, 
which are likely to be sampling dust emission from dif¬ 
ferent environments within each pixel. 

These caveats notwithstanding, the spinning PAH 
model for the AME predicts: 

1 . A linear correlation of the PAH surface density and 
the AME flux density at 30 GHz. 

2 . Relatively constant AME per PAH surface density 
over a range of radiation field strengths. 


3.2. Data Model 


T 353 is equal to the product of the dust mass column 
density and the dust opacity at 353 GHz. Because dust 
grains are much smaller than the wavelength of light at 
this frequency, the dust opacity is insensitive to the size 
distribution of the grains and 7353 is thus a reliable in¬ 
dicator of the total dust column density. There is evi¬ 
dence, however, that the dust optical pr operties change 


somewhat in diff erent environments (|Planck Collabora¬ 
tion et ah] |2014d). 


'The dust radiance 7 Z is the frequency-integrated dust 
intensity. TZ is estimated using the best-fit single temper¬ 
ature modified blackbody with r v oc z/^, and is a tracer 
of both the dust column density and the strength of the 
radiation field heating the dust. 

Finally, the 12 /am flux density is effectively a measure 
of the total power emitted by PAHs, i.e., the starlight 
power absorbed by PAHs. Thus, the ratio of the 12 /im 
emission to the dust radiance is a proxy for the fraction 
of dust in PAHs (see Equation^. The product /pah ^353 
is then a proxy for the PAH column density. It is this 
quantity which the spinning dust model predicts to be 
the most accurate predictor of the strength of the AME. 

We assume that t 353 , /paht 353 , i’ 2,ln \ and 


correlate with the AME intensity in a linear way, i.e., for 
each pixel i and each observable Ai , the AME intensity in 
that pixel is given simply by where cq is a constant 
to be determined. 

To identify the physical quantity which is the best pre¬ 
dictor of the AME intensity, we wish to quantify the in¬ 
trinsic dispersion in cg across all pixels. We assume that 
each pixel samples a Gaussian distribution with mean a 
and standard deviation a a . 

The likelihood of this model given the data over all 
pixels is: 


Coc U 


\f 'T + 


: xexp 


/rAME 

30 GHz,i 


— aAi 


\2n 


0/i + (Ai<j a ) 


( 3 ) 

where 07 ,; is the uncertainty on U^OGHz r We thus 
seek the maximum likelihood values of a and a a for each 
observable A. In practice, we employ t he emcee Markov 
chain Monte Carlo code (Foreman-Mackey et al. 2013) 
to derive the best fit values and confidence intervals for 
each parameter assuming uninformative priors. 

If there are zero point errors in any of the maps used 
in this analysis, then GHz = aA + b is a more ap¬ 
propriate functional form for the relationship, where b is 
a parameter to be fit. We find that the introduction of 
an intercept has little impact on the results, such as the 
ordering of the cr a /cq or the value of a. 


4. RESULTS 

4.1. Correlation with PAH Abundance 

We test the predicti ons o f the spinning PAH hypoth¬ 
esis laid out in Section |3T] by relating the observational 
data to physical pro perti es of the dust through the model 
described in Section lT 2 l 

In Figure [2] we plot the AME flux density at 30 GHz 
against 7353 , /pah^ 353 , the dust radiance 7£, and the 
100 /im flux density /I 00 / 6111 . All four correlate highly 
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Fig. 2.— The 30 GHz AME intensity is plotted against (a) T 353 , (b) /PAHT 353 , (c) 7 Z, and (d) l / 00 . We divide the plot area into 
hexagons of equal area in log-space and color each according to the number of pixels that fall within that hexagon. Each panel has the 
same logarithmic area allowing for straightforward comparisons between panels. In doing this, we have restricted the range of each plot to 
exclude some outlying points. In each panel we also plot (solid blue) the line with slope equal to the best-fit value of a. All four quantities 
are excellent tracers of the AME, but it is clear that 7 Z traces the AME with the greatest fidelity and least dispersion. 


TABLE 1 

Correlation Analysis 


A 

Ol =b a a. 

CTc/ol 

a A /A (0.2) 


8240±1210 MJysr- 1 

W m~ z sr ~ 1 J 

0.19 

0.01 

r 12 fj, m 

(7.48 ± 1.66) x 10“ 3 

0.22 

0.09 

GOO /Ltm 

J-u 

(2.88 ±0.65) x 10“ 4 

0.23 

0.10 

T353 

219 ± 72MJysr -1 

0.33 

0.27 

/PAH^353 

1123 ± 487MJysr -1 

0.43 

0.39 


Note. — a and <j a are the best fit values ob¬ 
tained from Equation [3] for each observable A such that 
GHz = a A- The formal uncertainties on a and a a 
are less than the quoted accuracy in the table, generally 
of order 0.1%. cr^/A (0.2) denotes the value of <ja/A 
such that <J a /a = 0.2, indicating the level of observa¬ 
tional uncertainty needed in A to produce a relationship 
with /AME gH z comparably tight as that observed with 

n. 

with AME as expected, and we present the fit slope of 
the relation with each in Table [l] The tightness of each 
correlation is indicated by oWa', with 7 Z having the tight¬ 
est correlation (see Figure ^p). 

If the AME arises from spinning PAHs, we would ex¬ 
pect the emission to correlate better with /pah ^353 than 


t353. While it is clear that both are good tracers of AME, 
/pahi~353 has a larger dispersion about the best-fit rela¬ 
tion than 7353. Maximizing Equation [3j A = 7353 yields 
a relation with cr a /a = 0.33 while A = /pahi ~353 yields 
cr a /a = 0.43 (see Table fTI) . Thus, /pah does not appear 
to contain additional information about the strength of 
the AME not already present in 7353. 

The spinning PAH model predicts that variations in 
the AME intensity per unit dust mass should arise from 
variations in the abundance of small grains. Therefore, 
as a second test of the link between the AME and PAHs, 
we look for correlations between /pah and the AME per 

^353- 

In Figure |3^i we plot the AME intensity normalized by 
ar 3 53 against /pah, but we find no evidence for the ex¬ 
pected positive correlation. We quantify this with the 
Spearman rank correlation coefficient r s which, unlike 
the Pearson correlation coefficient, does not assume a 
functional form for the relationship between the two vari¬ 
ables. We find r s = —0.15, suggesting anti-correlation. 

Similarly, quantifying the correlation between /pah 
and the AME intensity normalized instead by a7Z (see 
Figure [ 3 ) 3 ) yields r s = —0.02, suggesting that /pah car¬ 
ries no information on the AME intensity not already 
present in the dust radiance. 

We note that /pah is itself correlated with both 7353 
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Fig. 3.— /pah is plotted against the AME intensity normalized by (a) CHT 353 and (b) alZ. Here and in subsequent plots, the isodensity 
contours corresponding to 25, 50, and 75% of the pixels enclosed are plotted in green. There is no apparent correlation between /pah and 
the AME, at variance with the spinning PAH hypothesis. 


and 1Z as demonstrated in Figure Hi /pah (measured on 
~ 1° scales) varies from ^0.15 —"0.30 over most of the 
unmasked sky, so that variations in the PAH abundance 
might have been expected to account for a factor of ^ 1.5 
dispersion in AME intensity per dust mass. It is therefore 
striking that we detect no correlation of AME/77, with 
/pah- 

The observed strong correlation between /pah and r353 
is evidence for PAH destruction in the diffuse ISM (low 
7353 )- Depletion of PAHs by a factor of ~ 3 to sup¬ 
press photoelectric heating in the warm ionized medium 
(WIM) has been invoked t o explain high ratios of H a to 
free-free emission (Dong & Draine||2011) and is roughly 
consistent with the range of ,/pah we observe. _ 

A study of the Ha-correlated AME by |Dobler et al.| 
(|2009|) found the AME to be a factor of ^ 3 less in¬ 
tense than th e amplitude calculated by |Draine fc Lazar- 1 


than t 353. We would expect 7 Z and 7353 to be related 
through the strength of the radiation field- a fixed quan¬ 
tity of dust will radiate more when exposed to more ra¬ 
diation. The most straightforward conclusion is that the 
AME is enhanced by a stronger radiation field, which 
runs counter t_o t he p redictions of the spinning dust hy¬ 
pothesis (Ah-Haimoud et al. 2009 Ysard & Verstraete 


(1998a|b) and suggested that this discrepancy was 
due to PAH depletion in the WIM. However, subsequent 
studies have suggested an AME amplitude a factor of 
^ 3 lower than the earlier estimates even in Galactic 


clouds (Tibbs et al. 2010 2011 Planck Collaboration 


2010 ). 

A positive correlation between the radiation field 
strength and the AME intens ity has been noted in both 
the Perseus molecular cloud (Tibbs et al. | 201 1 ) and the 
HII region RCW1 75 (|Tibbs e t al. 2012b) . This trend 
was also noted by Planck Collaboration et al. (2014b), 
who attribute the correlation to a positive correlation 
between the radiation field and the local gas density. 

The relationship between the AME/77, and XJ, illus¬ 
trated in IT> is best-fit by a power-law of index -0.97. 
If the AME were thermal emission in the Rayleigh-Jeans 
limit, we would expect I^ ME to scale linearly with X^, i.e 

^Ioghz « r 30 T d . Since 7^ oc X d +/3 if r v oc v? , I^ ME /1Z 
should therefore scale as T^~ 3 


h—4.65 • 


if the AME is 


dtl1|2,]l4b|T Thus, the reWte ot |DobIer ct aI.| ( l200'4 thermal !™ s 5 sion and 3 ~ L65 ; Fittin S the data with a 


indicate instead a robustness of the AlVIE strength across 
environments and, together with our results, further ev¬ 
idence that the AME is not associated with PAHs. 

The correlation between /pah and 7353 also sheds light 
on the apparent negative correlation observed between 
/pah and I^ ME /rs 53 in Figure [ 3 ^ 1 . Figures [ 2^1 and [5^ 
demonstrate that <^7353 tends to underpredict the true 
AME intensity at low values of 7353 and to overpredict at 
high values. Since /pah is positively correlated with 7353, 
it is not surprising that assuming a linear relationship 
between the AME and /pah^"353 only exacerbates those 
discrepancies. 

4.2. Correlation with the Radiation Field 

Figure [ 2 ] and Table [l] also indicate the surprising re¬ 
sult that the AME is more tightly correlated with 7 Z 


power X/ ‘ yields a substantially worse y^ ed ( 10.2 vs. 
4.6). The data suggest then that the AME is nonthermal 
emission, unless the dust opacity at 30 GHz has a steep 
temperature dependence (e.g. 730/7353 oc Xj ,T , where 730 
is the dust optical depth at 30GHz.). We caution that 
systematic effects arising from fitting a single -Td mod¬ 
ified blackbody to the thermal dust emission may also 
alter the correlation with XJ and we thus cannot rule out 
a thermal emission mechanism completely based upon 
these data alone. 

If the fit AME component is contaminated with emis¬ 
sion from other low-frequency foregrounds, then this may 
also induce correlations with Td. Indeed, Figure [ 6 ] re¬ 
veals some coherent large-scale structures in the map of 
AME/77, that are likely relate d to stro ng synchrotron 


emission (see Planck Collaboration et al. 2015b Fig- 
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Fig. 4. — Plotting /pah against both (a) 7-353 and (b) 7 Z, it is clear that a strong positive correlation is present with both, consistent with 
depletion of PAHs in the diffuse ISM. The majority of the pixels have 0.15 < /pah < 0.30, suggesting that fluctuations in PAH abundance 
could account for a scatter of a factor of at most ~ 1.5 in AME intensity unit per dust mass. 
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Fig. 5.— The 30 GHz AME intensity norm alized by alZ is plotted against (a) Td and (b) (3 determined by modifie d blackbody fits to the 
dust SED (|Planck Collaboration et al.|2014a|), and against T 353 determined from component separation (Planck Collaboration et al. 2015b) 
in panel (c). Since FAHs are depleted m dense regions, the spinning PAH hypothesis predicts that the AME per '/<, sEould be smaller m 
denser regions. These regions are likely to have more cold dust and thus a smaller /3, but no correlation is observed between (3 or 7353 
and the AME intensity per 1Z. The data do suggest a possible correlation with T^. We find a best-fit power-law of T^ 0 ' 97 (blue solid) 
and plot also the best fit relation assuming the AME is thermal emission with T 30 oc 7353 , i.e. /^ ME oc T 353 T^, hence I^ ME /R oc T ^ -4 ' 65 
(black dashed). The latter relationship is disfavored relative to the form er ( x^ H = 1 0.2 vs. 4.6) . Finally, we note that our conclusions are 
unchanged when we redo this analysis using T ^ and (3 determined by [Planc k Collaboration et al.| (2015b). 
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Fig. 6 .— The Mollweide projection full-sky map of the 30 GHz AME intensity normalized by aTZ. Although some large-scale features 
associated with strong synchrotron emission are present, overall the map has little correlation with the synchrotron intensity (r s = 0.11). 
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ure 46) and perhaps the Fermi bubbles (S u et _al . 2010). 
We discuss this possibility further in Section [5j though 
aside from a few clear structures, contamination appears 
to be minimal. 

For pixels containing regions with both high and low 
radiation intensities, fitting the A > 100 /im emission by 
a single-temperature modified blackbody will lead to sys¬ 
tematic errors, tending to overestimate T^ and underes¬ 
timate both /? and T353. Thus, if the correlation between 
the AME and the dust radiance is driven by the corre¬ 
lation between the radiation field intensity and the local 
gas density, we might expect the AME per 7Z to correlate 
with /?. We find no evidence for such a correlation (see 
Figure [SJ)). 


4.3. Correlation with 

The ratio of /^qghz t° /xm is often quoted in the 
literature. We obtain a value of (2.88 d= 0.65) x 10 -4 
(see Table [l]), consistent with other determinations. For 
instance, performing component separation on WMAP 
observations at intermediate Galactic latit udes, jDavies 


et al. 


Alves et al. 


(2006) derived a ratio of 3 x 10 4 
(2010) likewise find a ratio of 3 x 10 -4 combining obser¬ 


vations of radio recombination lines in the Galactic plane 
with WMAP data. Studying HI I regions in the Galactic 
plane with the Very Small Array at 33 GHz, |TodorovIc| 
(2010) find a ratio of 1 x 10 -4 . Analyzing a sample 
Gal actic clouds with Planck observations, |Planck| 


et al. 
of 


Collaboration et al. (2014b) derived a ratio of 2. 5 x 10 4 
Thu s, the AME com ponent identified by Planck Collabo 


ration et al. (2015b) has a strength relative to the 100 /im 


dust emission in good agreement with what has been ob¬ 
served in other studies. 

However, as discussed in detail by Tibbs et al. (2012a), 
this ratio is subject to significant variability. While the 
100 /im emission is a reasonable proxy for the total dust 
luminosity, its sensitivity to the dust temperature intro¬ 
duces significant non-linearities in the relationship. This 
effect is evident at particularly low and high values of 
jioo/xm pig Ure |2ji as evidenced by the somewhat non¬ 
linear shape of the scatter plot. As expected, the tight¬ 
ness of the correlation between the AME and 100 /am 
increases when the analysis is restricted to a set of pixels 
with similar X^. 


4.4. Correlation with /i 2 ^ m 

Finally, we find that the 12 /am emission is also tightly 
correlated with the AME with dispersion only slightly 
less than that of 7 Z. While it is tempting to read this as 
a vindication of the spinning PAH model, the foregoing 
analysis suggests that this tight correlation is merely the 
product of the 12 /am emission being an excellent tracer 
of both the dust column and the radiation field strength. 
/pah, he. the 12 /im emission per unit 7£, does not cor¬ 
relate with the AME intensity (see Figures^ and [5]). 

Thus, while the PAH emission is an excellent predictor 
of the AME strength, this appears to be by virtue of 
being an excellent predictor of the dust radiance rather 
than the result of an inherent link between the AME and 
PAHs. 


4.5. Emission from Magnetic Dust 
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Fig. 7. — As in each panel of Figure [ 2 ] we plot the WISE 12 p,m 
intensity against the 30 GHz AME intensity. The correlation is 
comparably tight as observed with 7 Z. 

If the AME is not spinning PAH s, co uld it be 


sion from magnetic grains? Draine & Lazarian (1999) 
predicted that, unlike spinning PAH emission, magnetic 
dust emission would be equally strong per dust mass in 
both dense and diffuse regions. Since the dust in dense 
regions will be cooler than that in diffuse regions, pix¬ 
els with significant dust emission from both diffuse and 
dense regions will have broader SEDs and thus are fitted 
by smaller values of (3 relative to diffuse regions when 
fitting the SED with a modified single T^ blackbody. 
Thus, the AME per dust mass is predicted to be neg¬ 
atively correlated with f3 in the spinning PAH model and 
uncorrelated in the magnetic dust model. In Figure^ 
we demonstrate that AME/77, and (3 are largely uncorre¬ 
lated (r s = 0.03). However, magneti c grains should emit 
thermally ( Draine fe Hensley]|2013 ), whereas the rela¬ 
tionship derived above between Ip Mtj /7Z and T^ favors a 
non-thermal emission mechanism. In addition, thermal 
emission from magnetic d ust seems likely to be polarized 
(Draine & Hensley||2013|) whereas observations find the 
AME to be minimally polarized (see Section [ 5 ] for further 
discussion of polarization observations). 


4.6. Dependence on Masks 

The lack of correlation between /pah and the AME in¬ 
tensity per 7353 or 1Z is a potentially serious problem for 
the spinning PAH hypothesis. We thus test the sensitiv¬ 
ity of this result to the region of the Galaxy examined. 
In Figure [51 we perform the same analysis as in Figure^ 
and quantify the degree of correlation with the Spearman 
correlation coefficient r s . 

Starting with the N = 51579 pixe ls re maining un¬ 
masked following the cuts discussed in |2.5[ in Figure [ 8^1 
we consider only the pixels with Galactic latitude \b\ > 
15°. The behavior very much mimics that observed in 
Figure [ 3 ) 3 , with no compelling evidence of a correlation. 
In Figure [5b, we examine pixels with \b\ < 15°. These 
pixels have more PAH emission per dust radiance (i.e. 
higher values of /pah) than those at higher latitudes, 
as would be expected from PAH destruction in diffuse 
regions. Again, however, there is no indication of a cor¬ 
relation of the AME intensity per 7 Z with /pah- 

We next examine in Figure [ 8)2 only those pixels in 
which one or both of the AME components are signif- 
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Fig. 8. — As in Figure [3] /pah is plotted against the 30 GHz AME intensity normalized by oflZ. We examine four different masks: (a) 
considering only higher Galactic latitudes |6| > 15°, (b) considering only lower Galactic latitudes |6| < 15°, (c) considering only pixels in 
which at least one AME component was significant at > 5cr, and (d) using the standard mask only. For reference, the dashed lines mark 
the median values in the standard mask. In all cases, there is no evidence for positive correlation between the AME and /pah- 


icant at greater than 5cr. The cut on AME significance 
does not change significantly the behavior observed in 
Figure [ 8^1 other than to eliminate some pixels with low 
/pah values. As illustrated in Figure [4j the pixels with 
low /pah also tend to have low surface brightness, which 
may be responsible for inducing the slight negative corre¬ 
lation observed. Finally, Figure | 8 ]i shows for comparison 
the same analysis performed on the standard mask (Fig¬ 
ure [3^). 

The non-correlation of the AME intensity and /pah 
is therefore robust to assumptions either on the AME 
significance or the region of the sky analyzed. We now 
turn to the implications of this result in the following 
section. 


5. DISCUSSION 

One of the largest sources of uncertainty in this analy¬ 
sis is the fidelity of the AME spectrum recovered from de¬ 
composition. Due to a lack of data between the WMAP 
23 GHz band and the Haslam data at 408 MHz, it is dif¬ 
ficult to constrain the relative contributions of the AME, 
synchrotron, an d free-free in the frequency range of in¬ 
terest. Planck Collaboration et al. (2015b) notes that 


their estimates of synchro tron em ission are lower than 


the 9yr. WMAP analysis (Bennett et al. 2 013 ] ) by about 
70% at high Galactic latitudes and factors of several in 
the Galactic plane, with the AME component estimate 


being larger as a result. It is therefore possible that the 
AME map has a non-negligible synchrotron component. 

The presence of contamination in the AME map, such 
as free-free or synchrotron, would of course mean that 
/pah is not related to the inferred AME intensity in a 
perfectly linear way even if the AME comes from spin¬ 
ning PAHs. Further, since /pah correlates with r 353 and 
7 Z, it is also expected to correlate with the free-free and 
synchrotron emission, which increase with increasing gas 
(and therefore dust) column density. 

As illustrated in Figure [%, I^ ME /alZ is weakly anti¬ 
correlated with the 30 GHz Tree-free emission per 1Z ( r s = 
-0.07). 

In Figure IT we plot I^ ME /a1Z against the 30 GHz 
synchrotron intensity per TZ. The two quantities show 
evidence of a weak positive correlation ( r s = 0 . 11 ). 

In Figure IT we plot I^ ME /oflZ against the CO(l-O) 
line intensity per 7 Z. The two quantities show possible 
evidence of a weak positive correlation ( r s = 0.05). 

It is likely that some degree of contamination from 
other low-frequency components exists in the AME 
maps. This is to be expected given the theoretical diffi¬ 
culties modeling their spectra and the relative lack of ob¬ 
servational constraints at very low frequencies. In the fu¬ 
ture, more detailed decomposition enabled by additional 
data will improve the fidelity of the inferred AME spec- 
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Fig. 9.— The 30 GHz AME intensity normalized by aTZ is plotted against (a) the 30 GHz free-free intensity, (b) the 30 GHz synchrotron 
intensity, and (c) the CO(l-O) line emission all normalized by 1Z. There is evidence for weak correlation in all panels. 
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Fig. 10.— As in Figure[8] /pah is plotted against the 30 GHz AME intensity normalized by olIZ. We examine four different sets of pixels: 
(a) pixels with T& between 19 and 21 K (typical of the high-latitude sky), (b) pixels with between 1.55 and 1.65, (c) pixels meeting 
both the criteria of panels (a) and (b), and (d) using the standard mask only. For reference, the dashed lines mark the median values in 
the standard mask. Thus even when the conditions of the local environment are held relatively fixed, there is no evidence for a positive 
correlation between the AME and /pah- 
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trum and test the conclusions derived here. We note that 
the observed correlations between the AME/7 Z and the 
other low frequency components are rather weak, mak¬ 
ing it unlikely that these components are driving the ob¬ 
served non-correlation between the AME and /pah- 

Even assuming no contamination in the fit AME com¬ 
ponent, correlations with the synchrotron emission are 
plausible. For instance, if the AME arises from spinning 
ultrasmall grains, it might be affected in synchrotron- 
bright supernova remnants where it may be enhanced 
if shattering in grain-grain collisions increases the ultra¬ 
small grain population, or suppressed if ultrasmall grains 
are destroyed by sputtering. 

The theoretical uncertainties in the models for these 
emission components underscore the need for obtaining 
ancillary data at lower frequency. Upcoming 5 GHz ob¬ 
servations from the C-Band All-Sky Survey (C-BASS) 
(King et al.||2014) and 2.3 GHz observatio ns from the S- 
Band Polarized All-Sky Survey (S-PASS) (|Carretti et al.| 
z009) will play an invaluable role in disentangling the 
low-frequency components. 

A second source of uncertainty is the ability of /pah = 
A (z/7/) 12/xm /A7 Z to trace the PAH abundance. In par¬ 
ticular, the fraction of the PAH emission appearing in the 
WISE 12 /i m band can depend on the ionization state of 
the PAHs and other properties of the local environment 
(e.g. Draine & Li 2007 Draine 2011). Thus it may be 
possible to “wash out” a correlation between /pah and 
the AME intensity even if the AME arises from spinning 
PAHs. 

To test the plausibility of such a scenario, in Figure [l0| 
we correlate the AME intensity per 1Z with /pah only 
in pixels with similar environmental conditions as de¬ 
termined by their fit dust temperature and /3. We find 
no significant differences from our primary analysis and 
thus no evidence that variations in the local environment 
are driving the lack of correlation between /pah and the 
AME intensity per 7 Z. 

In addition to the uncertainties discussed above, each 
observable A also has an associated observational uncer¬ 
tainty a a which we have not included in Equation [3] a a 
is highly degenerate with cr a and thus it is difficult to 
quantify their relative contributions to the total uncer¬ 
tainty. Instead, we ask what assumed fractional uncer¬ 
tainty on A would be required for the data to be consis¬ 
tent with an intrinsic relationship with a a /a = 0.2, com¬ 
parably tight as found with 7 Z. We denote this quantity 
cja/A (0.2) in Table [11 finding that for t 353 and /paht 35 3 , 
cr a/ A would need tooe larger than 25% for the observa¬ 
tions to be consistent with an intrinsic relationship with 
so narrow a dispersion. The small uncertainties on the 
thermal dust emissio n at 353GHz reported by (Planck 
Collaboration et al.| 2015b) are inconsistent with what 
would be required Tor the data to be compatible with a 
tight intrinsic relationship between AME and either r 3 5 3 
or /pah735 3 - Thus, it does not seem plausible to attribute 
the entirety of our findings to the relative uncertainties 
of the observables. 


and 7 Z map s is inappropriate for the reasons discussed in 
Se ction|2.3[ it remains possi ble thatthe 7 Z map produced 
by |Planck Collaboration et al.| (|2014a|) is biased in ways 
that affect the correlation with /pah- Dust modeling 
employing additional data, such as 60 /im dust emission, 
will help clarify this issue. In the meantime, we note 
that 7£, being the integral of the FIR dust intensity, is 
rather insensitive to the specifics of modeling as long as 
the model provides a reasonable fit to the data. Thus, 
particularly when 100 /im data is employed, 7 Z is straight¬ 
forward to determine and we have no reason to suspect 
strong systematic biases. 

If indeed the AME is not correlated with the PAH 
abundance and is correlated with the strength of the ra¬ 
diation field, what are the implications for the origin of 
the emission? We present two possibilities: 

1 . The AME is spinning dust emission that arises pri- 
m arily from u ltrasm all grains that are not PAHs. 


Li & Draine (2001) have shown that as much as 
10 % of the Interstellar silicate mass could be in 
ultrasmall grains. However, it remains to be seen 
whether these grains could produce enough rota¬ 
tional emission to account for the AME. 


2 . The AME is predominantly thermal dust emis¬ 
sion, such as magnetic dipole emission from mag¬ 
netic materials. However, the relationship be¬ 
tween AME/TZ and T^ does not appear consistent 
with thermal emission (see Figure [5]). Addition¬ 
ally, current models of magnetic dipole emission 
do not predict behavior that could emulate the ob¬ 
served AME SED without invoking highly elon¬ 
gated (e.g., 5:1 pr olate spheroids) Fe inclusions 
(Draine & Hensley 2013|). Furthermore, thermal 
dust emission from an aligned component of the 
grain population is expected to be significantly po¬ 
larized, including magnetic di pole emissio n from 
ferromagnetic inclusions (Draine & Hensley |20 13), 
at odds with current non-detections of AME po- 
larization. Nevertheless, the theoretical calcula¬ 
tions of the emissivities of these magnetic mate¬ 
rials are still quite uncertain, and more laboratory 
data is needed to assess the behavior of magnetic 
materials at microwave frequencies to determine 
whether such grains could be a potential source of 
the AME. It remains conceivable that some inter¬ 
stellar grain component might produce the AME 
by non-rotational electric dipole radiation, but we 
are not aware of materials that could do this. 


Finally, Planck Collaboration et al. (2015c) demon¬ 
strate that the correlation between AME/ 1Z and /pah 
can be affected by systematic errors in the determination 
of 1Z. Although an analysis using the Commander-derived 
7 Z and /pah determined by straight division of the W?> 


Invoking an alternative explanation for the AME also 
requires explaining why the PAHs are not a substantial 
source of spinning dust emission. Because the PAHs are 
clearly present and must be rotating, this would require 
that the electric dipole moments of the PAHs have been 
significantly overestimated (the spinning dust emission 
scales as the square of the dipole moment). The elec¬ 
tric dipole moments of selected hydro carbon molecules 
compiled by Draine & Lazarian (1998b) have a scatter of 
nearly an order of magnitude. Further, harsh UV irradia¬ 
tion more easily destroys asymmetric molecules, perhaps 
preferentially selecting for a population of more symmet¬ 
ric PAHs with smaller dipole moments. Thus it is plau¬ 
sible that the dipole moment distributions adopted in 
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spinning dust models may significantly overestimate the 
tr ue electric dipole moments of interstellar PAHs. 


Lazarian & Draine (2000) estimated the polarization 
fraction of spinning dust emission to be p < 0.01 near 
30 GHz. This estimate would also apply to spinning dust 
emission from non-PAH grains. However, if the AME is 
thermal emission from large aligned grains, it should be 
significantly polarized with E _L B 0 for electric dipole 
radiation or E || B 0 for magnetic dipole rad iation from 
magnetic inclusions ( Draine fe Hensley||2013 ). 

Observations of known A 1 V 1 E sources m polarization 
suggest minimal polarization of the AME. In the Perseus 
molecular cloud, the polarization fraction of the total 


emission was found to be at 11 GHz (Battistelli 

|2006h , less than 6.3% at 12 GHz (|Genova-Santos 



2015;), less than 2.8% at 18 GHz (Genova -bantos 
|2015p, an d les s than ~ 1% at 23 G Hz 


3 GHz (|L 6 pez- 

)1 lb. Likewise, 


Garaballo et al. 2011 Dickinson et al. 2011). Likewise, 
observations of the p Ophiuchi molecu lar cloud have 
yielded upper limits of < 1 % at 30 GHz (|Casassns et al 


et al. |[20ll]). |Mason et al.|(|2009| placed 
: 2.7% on the 9.65 GHz polarization frac- 


2008 Dickinson 
an upper limit of 
tion in the dark cloud Lynds 1622. 21.5 GHz observations 
of the Hu regi on RCW175 yielded a polarization fraction 
of 2.2 ±0.4% (|Battistelli et al.|2015|) , though it is unclear 
whether this polarization is arising from the AME or a 
sub-dominant synchrotron component. 

These upper limits on polarization in the 10 - 30 GHz 
emission appear to favor spinning dust emission from a 
non-PAH population of ultrasmall grains. 


6 . CONCLUSION 

We have combined the Planck foreground component 
maps, Planck modified blackbody dust parameter maps, 
and WISE 12 /im maps to test key predictions of the 
spinning PAH hypothesis. The principal conclusions of 
this work are as follows: 


1 . t 353 , the dust radiance 7and /i 2 ^ m are all excel¬ 
lent predictors of the 30 GHz AME intensity. 7 Z ex¬ 
hibits the tightest correlation, suggesting that the 
AME is sensitive to the strength of the radiation 


field. 

2 . Neither AME /7353 nor AME/77. show any correla¬ 
tion with the PAH emission whether considering 
the full sky, regions close to the Galactic plane, or 
higher Galactic latitudes. 

3 . We find that /pah is correlated with both 7353 and 
7consistent with PAH destruction in low density 
regions. 

4. Taken together, these facts pose a serious chal¬ 
lenge to the spinning PAH paradigm as the expla¬ 
nation for the AME. Alternative explanations, such 
as magnetic dipole emission from ferro- or ferri- 
magnetic grains, should be more thoroughly inves¬ 
tigated. 

5. More low frequency constraints are needed to break 
degeneracies between the AME, free-free, and syn¬ 
chrotron to enable more accurate decomposition 
and to better constrain the AME spectrum. Up¬ 
coming all-sky observations from C-BASS and S- 
PASS will thus facilitate deeper investigations into 
the origin of the AME. 

6 . Further measurements of AME polarization will 
help clarify the nature of the grains responsible for 
the AME. 
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